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Abstract 

The governing equation is u t — (a(x)u x ) x , < x < 1, t > 0, u(x,0) = 0, u(0,t) = 
0, a(l)u'(l,t) = f(t). The extra data are u(l,t) = g(t). It is assumed that a(x) is a 
piecewise-constant function, and / ^ 0. It is proved that the function a(x) is uniquely 
defined by the above data. No restrictions on the number of discontinuity points of 
a{x) and on their locations are made. The number of discontinuity points is finite, 
but this number can be arbitrarily large. 

If a(x) £ C 2 [0, 1], then a uniqueness theorem has been established earlier for 
multidimensional problem, x G R n ,n > 1 (see MR1211417 (94e:35004)) for the sta- 
tionary problem with infinitely many boundary data. The novel point in this work is 
the treatment of the discontinuous piecewise-constant function a{x) and the proof of 
Property C for a pair of the operators {t\, £2}, where £j := + fc 2 g|(x), j = 1,2, 

and q"j(x) > are piecewise-constant functions, and for the pair {Li,L2}, where 
LjU :— — [a,j(x)u'(x)]' + Xu, j = 1,2, and a,j{x) > are piecewise-constant functions. 
Property C stands for completeness of the set of products of solutions of homogeneous 
differential equations (see MR1759536 (2001f:34048)) 

Keywords: inverse problems, heat equation, Property C, piecewise- constant ther- 
mal conductivity. 
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1 Introduction 

Let 

u=(a(x)u'Y, 0<x<l, i>0, u':=— , u:=—, 1) 

u(x, 0) = 0, u(0, t) = 0, a(l)u' (1, t) = f(t) £ 0, (2) 

u(l,t)=g(t). (3) 

Problem (JTJ) — (fZJ) describes the heat transfer in a rod, a(x) is the heat conductivity, a(l)n'(l, t) 
is the heat flux, g{t) is the measurement, the extra data. 
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The inverse problem (IP) is: 

IP: Given f(t) and g{t) for all t > 0, find a(x). 

Assumption A: a{x) is a piecewise- constant function, a(x) = aj, Xj < x < 
X\ = 0, x n+ i = 1, < Cq < dj < c\, 1 < j < n. 

This assumption holds throughout the paper and is not repeated. The set of piecewise- 
constant functions with finitely many discontinuity points is denoted by II. 

If a(x) £ C 2 , then the uniqueness of the solution to some multidimensional inverse 
problems has been proved in [3] (see also [2]). Problem {I])-© with a(x) € C 2 ([0, 1]) 
has been studied in py, [6]. The treatment of discontinuous piecewise-constant a(x) is of 
interest in applications. 

In p] equation dU) with the conditions u(0, t) = u(l,t) = 0, u(x, 0) = g(x), was studied, 
and the measured (extra) data were the values u(£ m ,t), Vi > 0, 1 < m < M, < £ m < 1, 
where M = 3n, and n is the number of the discontinuity points of a(x). It was assumed in 
p] that minj \xj — Xj+\ \ is not too small. Under these assumptions the uniqueness theorem 
for the IP was proved in [I], and an algorithm for finding a(x) was proposed. The stability 
of this algorithm with respect to perturbations of the data was not studied in pQ. 

In our paper the extra data ([3]) consists of measurement, taken at one point, rather than 
at 3n points, and we impose no restrictions on minj \xj — xj + i\. Under these assumptions, 
which are much weaker than in pQ , we prove the uniqueness of the solution to IP. 

One of our main results is 

Theorem 1 The IP has at most one solution. 

Remark 1 The IP is ill-posed: small variations of the data {f(t),g(t)} in the C(0, oo)— norm 
may lead to large variations of the coefficient a(x), or may lead to a problem which has 
no solutions. We assumed that the data are known for all t > 0. If one assumes that 
f(t) = for t > T, where T > is an arbitrary fixed number, then the solution u(x,t) is 
an analytic function oft in the region t > T. Therefore the data {f(t),g(t)}, known in 
the interval [0,T + e), where e > is an arbitrary small fixed number, determine uniquely 
the data for all t > 0. Thus, if f(t) = for t > T, then the uniqueness theorem for the 
solution to IP remains valid if the data are known for t £ [0, T + e). 

Let us formulate IP in an equivalent form. 

Take the Laplace transform of the equation dU)-©, denote 




and get: 



Xv 



(a(x)v')' = 0, 
a(l)v'(l,X) = 



0<x<l, u(0,A)=0, 
F(X), v(l,X)=G(X), 



(4) 
(5) 



where F := Lf and G := Lg. 

The IP can be reformulated as follows: 
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IP: Given F(X) and G(X) for all A > 0, find a(x). 

Let us transform equation ([H)-© to yet another equivalent form. 

Let a(x)v' := ip. Then (JH)-© can be replaced by the following problem 

-< + Aa- 1 (x)^ = 0, 4>(1,X)=F(\), ij/(0,X) = 0, (6) 

^(1,A)=AG(A). (7) 

The IP can be reformulated as follows: 

IP Given G(X) and F(X), find a -1 (a:) := q 2 (x). 

Let 

£tP ■- + k 2 q 2 {x)'ip = 0, X:=k 2 , q 2 {x) : = a _1 (x), cf/ 1 < g 2 (x) < c \ (8) 
Consider the following problems: 

£jiPj = 0, e j: =-^ + k 2 q 2 (x), tf j (0,k)=0, Vi(0,fc) = l, j = l,2. (9) 

Our second main result is 

Theorem 2 The sets {ipi(x, k)tp2(x, fc)}vfe>Q and A)i4(x, A)}va>o> & := A 1 / 2 , are 

dense in the set H of piecewise- constant functions on [0, 1]. 

Remark 2 Theorem^ says that if h(x) £ II and 

/ h(x)ipi(x,k)ip 2 (x,k)dx = 0, Vfc > 0, (10) 

i/ien /i = 0. Similar conclusion holds iftpj(x,k) is replaced by v'j(x,X) in ([10]) , 5W/i a 
property of the pair of the operators {£±,£2} is called Property C (f^j, J^]/j. 

Clearly if the set {ipi(x, k)ijj2(x, fc)}vfc>o is dense in the set II, then the set of products 
{v'i(x, X)v' 2 (x, A)}va>o * s dense in the set II. 

In Section 2 proofs are given. 

2 Proofs 

2.1 Proof of Theorem [1] 

Proof. We prove this Theorem for the problem dU)-©. Suppose there are Vj and aj £ II, 
j = 1,2, which solve problem (JU)-©, and let w :=v\ — V2- Then 

Xw - (aiw'Y = (pv 2 )' , p := ai(x) - a 2 (x), (11) 
tu(0,A) = 0, «;(1,A)=0, a lV [(l, X) = a 2 v 2 (l, A). (12) 
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Multiply (fTTj) by v%, a solution to equation with a = 01, and integrate over [0, 1], and 
then by parts, to get 



p(x)v' 2 v[dx 





1 


l 


pv' 2 vi 


+ aiizA/i 


— aiit^ 











0. 



VA > 0, X = k 2 , k> 0, (13) 



where we have used the conditions w(0, A) = w(l, A) = and ai (1)^(1, A) = a 2 (l)v' 2 (l, A). 
Note that V2(x,X) can be considered as an arbitrary solution to equation up to a 
constant factor. The set {v[(x, X)v' 2 (x, A)} is dense in II by Theorem [2] Since a\(x) — 
a,2 (x) := p(x) € II, it follows from (fT3l) that p(x) = 0. So a\ = ai- Theorem [T] is proved. 



□ 



2.2 Proof of Theorem H 

Proof. Let us prove completeness of the set of products (x, k)ip 2 (x, fc)}vfc>o- Assume 
that h £ II and (|10p holds. The function ipj(x,k), j = 1,2, are entire functions of k. This 
follows from the integral equation for ipj , which is an immediate consequence of equations 
©-©: 

ipj(x,k) = 1 + k 2 / (x - s)q](s)ipj(s,k)ds, x > 0, j = l,2. (14) 

JO 

Equation p^]) implies that for any fixed k one has ipi(x) := ipj(x,k) > 1, Vx £ [0,1], 

j = 1,2, that ^.(a?,*) > 0, t/j"(x,k) > 0, and &m ^' k) > for all m = 0,1,2, 

Consequently, ^j(x), j = 1,2, are convex functions of x on the semiaxis x > 0. Since 
ij)j(x,k), j = 1,2, are positive, it follows from (fT4j) that ipj(x,k), j = 1,2, are increasing 
functions with respect to both x and k. So we have 

ipj(x, k) > 0, ^(x, k) > 0, ip"(x, k) > 0, VA; > 0, j = 1, 2. (15) 

Assume < x\\ < x\i < • • • < xin- l < 1 and < x 2 i < x 2 2 < ■ ■ < x 2 n 2 < 1 are 
discontinuity points of a±(x) and a 2 (x), respectively. 

To derive from (fT0|) that h = it is sufficient to prove that h(x) = 0,Vx G [xo,l], 
where xo := max(xiAr x , x 2 n 2 )> because then one can prove similarly, in finitely many steps, 
that h = on the whole interval [0, 1] using the assumption h & II. We have 

ip'-(x, k) = k 2 q] Nj (x^j(x, k), Vk > 0, Vx G [x , 1], (16) 

where qjjf. is the value of qj on the interval [xq, 1]. From ([16]) one gets 

i/)j(x, k) = aj (k)e kq ^ {x - xo) + b j (k)e- kqiN ^ x - xo) , V/c > 0, j = 1, 2. (17) 
It follows from (TJS]) and ^ that 

^•(x , ft) = aj(k) + bj(k) > 0, ^(x , k) = kq jNj [ aj (k) - bj(k)\ > 0, (18) 
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and 



This implies 



2a,j(k) = i>j(xo, k) + 



kqjNj 



> i/jj(xo,k). 



aj(k) > \bj(k)\ > 0, VA;>0, j = 1,2. 
Since h E II, one may assume without loss of generality that 

h(x) = C>0, VxG[x ,l]. 

It follows from ([IOJ) that 

l>XQ 

Jo 

From (H5]), ([HD and ([20}, one gets 

1 < ipj(x, k) < ipj(x , k) < 2cij(k), < x < xq, VA; > 0, j = 1,2 
Therefore, 

f-XO 

<Aai{k)a 2 {k) \ \h(x)\dx. 



f 1 

ipi[x, k)ip2\X, k)h[x)dx = / ^i{x,k)i^2{x,k)h{x)dx, Mk 

J XO 



(19) 
(20) 

(21) 

> 0. (22) 



(x, k)ip2(x, k)h{x)dx 
From (|2"0"]). (fT7|) and {IS} one obtains 

ipj(x,k) > aj (k)[e kqjN ^ x - Xo) - e- kqjN ^ x - xo) ] 



x e [x ,l], j = 1,2. 



(23) 
(24) 

(25) 



Take an arbitrary y € (xq, 1) and fix it. One has iftj(x, k) > ipj(y, k), \/x € [y, 1]. Therefore, 



■01 (x, k)tp2(x, k)h(x)dx > C(l — y)ipi{y, k)ip2(y, Vfc > 0. 



■«'(> 



This, 



and (|24[) imply the following inequalities: 

Mv> k )Mv,k) 



oo > 4 / |/i(x)|dx > C(l-y)- 
Jo 



ai(k)a 2 (k) 



Mk > 0. 



It follows from (ESI) that 



lim — — — — = oo. 



(26) 



(27) 



(28) 



Let k — > oo in ([27]) and use ([28]) to conclude that C = and, therefore, /i(x) = for 
x € [xo, 1]. Similarly one proves that /i(x) = for all x £ [0, 1]. 

Theorem [2] is proved. □ 
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